COSMOLOGICAL PARAMETERS FROM COSMIC SHEAR 



MARTIN KILBINGER 

Institut fur Astrophysik, Universitdt Bonn 

Abstract. We present simulations of a cosmic shear survey and show how the survey geometry 
influences the accuracy of determination of cosmological parameters. We numerically calculate 

the full covariance matrices of the two-point statistics £+, £;_ and ^M| p \ and use maximum 
likelihood and Fisher matrix analyses in order to derive expected constraints on the cosmological 
parameters Cl m ,as,T,n s and J7a- 



1 Introduction 

For a thorough introduction to cosmic shear, see the review of Refregier [1] in this volume. 

Since its first detection in 2000, cosmic shear has proven to be a valuable tool for cosmology. 
Constraints on cosmological parameters, in particular the (dark+ luminous) matter density f2 m and 
the power spectrum normalization as has been obtained from cosmic shear surveys with observed 
of up to several dozen square degrees (see Table 1). 

A cosmic shear survey has to cover a large area containing hundreds of thousands of galaxies, whose 
shapes can be determined. Because telescope time is limited, one has to carefully choose the locations 
of the pointings, in other words, the geometry of the survey: On the one hand, many independent 
lines-of-sight are preferable, lowering the sampling or "cosmic variance". On the other hand, it is 
important to measure the shear on a large range of angular scales. Even with modern wide-field 
imaging cameras, separations of more than ~ 1 degree cannot be accessed by individual fields-of-view 
- one has to observe some fields near to each other and measure galaxy shape correlations across 
fields. Thus, a tradeoff between "clustering" and wide separation of pointings has to be found. In this 
contribution, we present a method to compare different survey configurations with respect to their 
ability to constrain cosmological parameters. For a more detailed description, see [2]. 

The survey geometries considered here consist of n images which are placed into P patches of 
radius R, each patch containing N pointings (thus n = P ■ N). The individual patches are assumed 
to be uncorrelated i.e. widely separated and randomly placed on the sky. The largest scale on which 
cosmic shear can be probed in these cases is 2R. In addition to these patch geometries, a survey is 
considered where all n individual fields are uncorrelated. For a given n, this survey has the smallest 
possible cosmic variance, but the largest scale to probe cosmic shear is only \/2 times the image size. 

We use a total image number of n = 300. The individual images are 13' x 13'-fields (which is 
roughly the field-of-view of VIMOS). This corresponds to a survey area of about 14 square degrees. 
The number of galaxies, for which a shape measurement is feasible, is set to 30 per square arc minute, 
corresponding to a limiting i?-band magnitude of about 25.5. 



Table 1: Recent and ongoing cosmic shear surveys. m\\ m is the limiting magnitude and zq is the 
redshift parameter of the source galaxy distribution, see also eq. (1). 
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Figure 1: The diagonal elements of 
the covariance matrices (2), split up 
into the individual terms D, M and 
V. Upper row: 300 single uncorre- 
cted images, where the largest scale 
is y/2 ■ 13' (thick lines) and a patch 
geometry with N = 60 and R = 80' 
(thin lines). Lower row: N = 10, 
R = 100' (thick lines) and N = 60, 
R = 100' (thin lines). 



1 10 

6 [arcminl 



1.1 1 10 100 

6 [arcminl 



2 Cosmological model 

For the power spectrum of the matter fluctuations, we assume an initial power law P\ oc k ns , the 
transfer function for Cold Dark Matter from [8] and the fitting formula for the non-linear evolution of 
[9]. The redshift distribution of the source galaxies is set to [10] 

where V denotes the Eulerian gamma function. Our reference cosmology is a flat ACDM model with 
= 0.3, f^A = 0.7, n s = 1, the shape parameter T = 0.21 and the normalization ag = 1. The 
parameters of the redshift distribution are zq = 1 and (3 = 1.5, which corresponds to a mean source 
redshift of ~ 1.5. 



3 Covariance of two-point statistics of cosmic shear 

We use three different second order statistics of cosmic shear for our analysis, the two two-point 
correlation functions £± and the aperture mass (m^J. 

The covariance matrices of the two-point correlation functions are defined as follows: 

Cov ±± = Cov(|±,0 i; | ± ,%) = ((i±(0i) - £±(04)) (L(0j) - £fc(0i))) , (2) 

where £± are unbiased estimators of £±. The covariance of the aperture mass is defined in complete 
analogy. 

The covariance matrices of the correlation functions consist of three different terms: a shot noise 
term (D) which appears only on the diagonal and depends solely on the intrinsic ellipticity dispersion, 
a term which depends on cosmic variance only (V) and a mixed term (M). In Fig. 1, the diagonal 
elements of the covariance matrices are plotted, split up into the individual terms. For Cov ++ , the 
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Figure 2: 1-cr, 2-a and 3-cr con- 
fidence contours of the figure-of- 
merit (3). The parameter Ss is pro- 
portional to as, see text. Solid lines 
are from the combination of and 
dotted (dashed) lines are calcu- 
lated using £ + (£_) alone. The left 
panel shows the results for the 300 
uncorrelated images, the right panel 
is for a (N = 60, R = 80')-patch ge- 
ometry. 



cosmic variance V dominates on scales larger than about one arc minute, whereas Cov is dominated 

by D. As shown in the upper row of Fig. 1, there is about a factor of two in V+ + between the two 
extreme geometries regarding cosmic variance. The lower row of panels of Fig. 1 compares patch 
geometries with the same radius, but with a low (N = 10) and a high (N = 60) image density in the 
patches. In the first case, there is a sharp transition at the image boundary scale, which is not present 
for the denser patches. 



4 Analysis 

4.1 Maximum likelihood 

Using the covariance matrices, we construct a figure-of- merit 

x 2 (p) = £ (6W " £) Cov ^ " 3) > 6 = (3) 

where the superscript * denotes the fiducial model (see Sect. 2) and p is a set of cosmological parameters 
which is tested against this model. This figure-of-merit can be calculated using the two correlation 
functions and £_ separately, or the combination of both; the resulting functions are called x\i x 2 - 
and x?ot> respectively. 

We calculated (3) for a number of different survey geometries. The most interesting constraints 
on cosmological parameters from cosmic shear are those on f2 m and a$. We kept all other parameters 
fixed and reparametrized to Eg = er 8 /[0.41 + 0.59 ($l m /0.3)~ ' 68 ] in order to compensate for the 
high elongation of the Q m -as- likelihood contours. 

In Fig. 2, the f2 m -E8-likelihood contours of (3) are plotted, for two extreme geometrical configura- 
tions regarding cosmic variance: the uncorrelated images and the (N = 60, R = 80')-patch geometry. 
Clearly, the x~ -contours are more extended in the case of the uncorrelated images than for the patch 
geometry. This is because £_ contains much information on large scales which is absent in the case 
of the uncorrelated images. In contrast, the x+- con ^ ours are tighter in this case than for the patch 
geometry. 

Furthermore, in both cases, the difference between x\ an d Xtot are small. Thus, most of the infor- 
mation concerning these two cosmological parameters is contained in £ + ; the additional information 
coming from ^_ is relatively small. 

4.2 Fisher information matrix 

The Fisher information matrix F is defined as the expectation value of the second derivative of the 
likelihood function. From the Cramer-Rao inequality, we get a minimum variance bound (MVB) 
a (pi) = V {F^^u for any unbiased estimator of a parameter pj. 




Figure 3: l-u-likelihood contours (solid lines) compared with the 1-cr-error ellipse from the Fisher 
matrix: the dashed ellipse is for a flat Universe (as it is the case for the likelihood contours), the 
dotted one is for Q\ = 0.7. a/b is the axis ratio of the ellipses (the case Q\ = 0.7 is in parentheses). 
The configuration is a survey of 300 uncorrelated images, we used the combined and 



In Fig. 3, we compare the MVB from the Fisher matrix with the 1-er-contours of the likelihood, 
using in both cases. As expected, the likelihood contours are larger than the 1-cr-ellipse from 
the Fisher matrix. The orientation of the Fisher ellipse coincides with the likelihood shapes, i.e. 
the direction of the minimal and maximal degeneracy of parameters is recovered. The larger the 
degeneracy between two parameters, the larger is the deviation between the local approximation by 
the Fisher matrix and the likelihood function. 

Next, we assume that trie parameters £^ m , cr 8 ,r and n s are to be determined simultaneously from 
the survey and calculate the MVB for these parameters, using different survey geometries. We now 
perform the analyses using the M ap -statistics which is thought to be the most useful statistics for 
cosmic shear because of its ability to separate E- from B-modes. 

The MVB for these parameters for a variety of survey geometries is plotted in Fig. 4. The smallest 
values for the MBVs are obtained for a configuration with N = 30 and R = 60' — 70'. Note the 
similarity of the curves for O m and as which both determine more or less the amplitude of the power 
spectrum, and those for T and n s which are responsible for the shape. 



5 Conclusions 



The best constraints on cosmological parametes are obtained for a survey with N = 30 images dis- 
tributed in patches on the sky. The dependence on the patch radius is quite small. The MVBs of 
some combinations of cosmological parameters to be determined from a cosmic shear survey with 
N = 30 and R = 60' are given in Table 2. The cosmological constant is only poorly constrained. For 
fi m , us, r and n s , the optimal survey consists of small patches of about R = 60 arc minutes of radius, 
with N = 30 images randomly distributed in each patch. The difference of the MVB makes up to 
25 percent for different survey geometries, thus a 25 percent improvement on the determination on 
cosmological parameters can be obtained solely by choosing the appropriate survey geometry. 
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Table 2: Each row shows the MVBs for a combination of three 
or four cosmological parameters which are to be determined 
simultanously from a cosmic shear survey. 
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Figure 4: The MVB for the parameters f2 m , T, as and n s , using the ^Mf p y-statistics. $7 m , T, as and n s 

are assumed to be determined from the data, all other parameters are kept fixed. Each point in the 
plot represents a patch geometry with N images in patches of radius R. 
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